Viscous flow over a chemically patterned surface 
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The classical fluid dynamics boundary condition of no-slip suggests that variation in the wetta- 
bility of a solid should not affect the flow of an adjacent liquid. However experiments and molecular 
dynamics simulations indicate that this is not the case. In this paper we show how flow over a 
solid substrate with variations of wettability can be described in a continuum framework using the 
interface formation theory developed earlier. Results demonstrate that a shear flow over a perfectly 
flat solid surface is disturbed by a change in its wettability, i.e. by a change in the chemistry of the 
solid substrate. The magnitude of the effect is shown to be proportional to cos Si — cos #2 where 
#1 and 62 are the equilibrium contact angles that a liquid-gas free surface would form with the two 
chemically different parts of the solid surface. 

PACS numbers: 68.08.-p, 47.61.-k, 68.43.-h 



I. INTRODUCTION 

The flow of liquids over chemically patterned surfaces 
is an exciting and relatively new area of fluid mechan- 
ics with applications in many emerging technologies 
Such flows are of particular interest in microfluidics where 
an increasing surface to volume ratio of liquids means 
that surface effects become of greater significance 
The correct description of the physics at solid-liquid in- 
terfaces then becomes imperative to the success of any 
attempt to model this class of flows It has been 
shown that by patterning a substrate with hydrophillic 
and hydrophobic regions it is possible to confine a liq- 
uid to a microchannel [J, [5( , to improve the accuracy of 
droplet deposition @, @] or to create a structured film § . 
Alternatively, when a wettability gradient is present, un- 
balanced surface tension forces can lead to the movement 
of liquid drops [1, [l(| • 

The effect of variable wettability of the solid substrate 
on the adjacent flow has been studied theoretically us- 
ing molecular dynamics simulations [ill [l2| . In this ap- 
proach, variations in wettability are modelled by varying 
an interaction potential between molecules of the solid 
and the fluid. The results show that a change in wet- 
tability does affect the flow, most notably producing a 
component of velocity normal to the solid surface. Given 
that "wettability" can be introduced as a macroscopic 
characteristic of a liquid-solid system, one should be able 
to model the effects discovered by molecular dynamics 
macroscopically using an appropriate formulation in the 
framework of continuum mechanics. 

In the continuum approximation, molecular proper- 
ties of the contacting media, as well as molecular length 
and time scales, do not appear explicitly. Instead, they 
manifest themselves in the bulk equations and bound- 
ary conditions via macroscopic transport coefficients and 
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parameters of constitutive equations. The classical no- 
slip boundary condition leaves no room for incorporating 
the effects of variable wettability and, as is well known 
[IH G3 > it becomes inadequate for dealing with processes 
of dynamic wetting where wettability of the solid plays a 
key role. 

A generalisation of no-slip often cited in the literature 
is the Navier slip condition [15|,ll6|, where slip, i.e. the dif- 
ference between the tangential velocities of the fluid and 
the solid, is assumed to be proportional to the tangential 
stress acting from the fluid on the liquid-solid interface. 
The coefficient of proportionality is referred to as the 'slip 
coefficient' or the 'coefficient of sliding friction'. Various 
expressions for this coefficient have been proposed [lTj 
some of which are based on the results of molecular dy- 
namics simulations [l8j . An appropriate choice of the 
slip coefficient allows one to describe situations in which 
the characteristic length scales of the flow are small (e.g. 
nanopores), so that the effects of slip become apparent 

The problem arises when the approach based on the 
Navier slip condition is used to describe flows over solid 
surfaces with variable wettability, since the existing mod- 
els provide no conceptual link between slip and wetta- 
bility. Indeed, as reviewed in §9 of [20| . in the area of 
dynamic wetting, where the physics of wettability be- 
comes a key factor, the so-called 'slip models' treat the 
behaviour of the contact angle, which is the macroscopic 
'measure' of wettability, and slip, which is used as a 
boundary condition removing the stress singularity at 
the moving contact line, as completely independent. On 
a macroscopic level, any coupling of a formula for the 
velocity-dependence of the dynamic contact angle and a 
slip boundary condition will produce a new slip model. 
Alternatively, linking the slip coefficient (or any other slip 
condition) to wettability via molecular dynamics simula- 
tions negates the advantages of the self-contained and 
experimentally verifiable continuum modelling. 

Applying the no-slip or constant slip coefficient Navier 
boundary condition at a solid surface implies that the 
flow of a liquid will not be affected by variations in the 
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wettability of an adjacent solid. This is in direct conflict 
with the conclusions of molecular dynamics simulations 
[IU HH that were specifically tailored to investigate this 
problem. In the present paper we consider flow over a 
solid surface of variable wettability using a model de- 
veloped earlier f2l| that incorporates wettability in the 
framework of continuum mechanics from the following 
perspective. Given that dynamic wetting is, by its very 
name, a process of creating a new/fresh liquid-solid in- 
terface ('a wetted solid surface'), an adequate mathemat- 
ical model of dynamic wetting must have the physics of 
interface formation process at its core. The model (here- 
after referred to as the interface formation model) has 
been derived from first principles [II[ using methods of 
irreversible thermodynamics and notably it ensures that 
all aspects of wettability (contact angle, slip, adsorp- 
tion) are interrelated. Embedding dynamic wetting in 
a general physical framework, as a particular case from 
a broad class of flows with forming/disappearing inter- 
faces, allows one to use the same model, without any ad 
hoc alterations, to describe other flows from this class 
[13, US IHl ■ This makes it possible to use the values for 
material constants determined from completely indepen- 
dent experiments. Comparison with such experiments 
have confirmed the model's validity and accuracy pH. I25I] . 

We will examine the fully representative case of a 
plane-parallel shear flow over a smooth solid surface that 
encounters a transition region between substrates of dif- 
ferent wettability. 

II. PROBLEM FORMULATION 

Consider the steady flow of an incompressible Newto- 
nian liquid passing over a stationary flat solid surface. 
The liquid is driven over the solid by a plane-parallel 
shear of magnitude S in the far field. We consider the 
bulk flow to be described by the Navier-Stokes equations 

V-u = 0, pu • Vu = -Wp + pV 2 u (1) 

where u, p, p and p are the fluid's velocity, pressure, den- 
sity and viscosity, respectively. The boundary conditions 
on a solid surface of uniform wettability, that follow from 
the interface formation theory f2~l| are given by 

pn-[Vu + (Vu)*] • (I - nn) + ±Vct = /?u • (I - nn), (2) 
pun = , (3) 

T 

V • (/v s ) = (4) 
r 

(v s - iu) • (I - nn) = aVa, (5) 
* = l(pto) - P S )- (6) 



Here a is the surface tension in the 'surface phase', i.e. 
a microscopic layer of liquid adjacent to the solid surface 
subject to intermolecular forces from two bulk phases; p s 
is the surface density (mass per unit area) of this layer 
and v s is the velocity with which it is transported; n is 
the unit vector normal to the solid surface pointing into 
the liquid; I is the metric tensor; (I — nn) is a tensor 
that extracts the tangential component of a vector, for 
example u • (I — nn) = um, where the double line sub- 
script denotes tangential component; a, /3, 7, r, p s , Q s and 
p s e are phenomenological material constants. 

The model has previously been discussed in detail (e.g. 
f2l|). so that here we will only briefly recapitulate the 
meaning of the terms. The surface tension is considered 
as a dynamic quantity related to the surface density via 
the equation of state in the 'surface phase' ^ which is 
taken here in the simplest linear form. The constant 7 
is associated with the inverse compressibility of the fluid 
whilst is the surface density corresponding to zero 
surface tension. As in Gibbs' theory of capillarity, the 
surface tension itself may be positive or negative depend- 
ing on whether the solid is hydrophobic or hydrophillic. 
Gradients in surface tension influence the flow, firstly, 
via the tangential stress boundary condition ([2]), i.e. via 
the Marangoni effect, and, secondly, in ([5|) by forcing the 
surface velocity v s , which is parallel to the solid surface, 
to deviate from that generated in the surface phase by 
the outer flow. The constants a and f3 characterise the 
response of the interface to surface tension gradients and 
an external torque, respectively; in the simplest variant 
of the theory both are properties of the fluid and have no 
relation to the wettability of the solid. Mass exchange 
between the bulk and surface phases, caused by the pos- 
sible deviation of the surface density from its equilibrium 
value pi, is accounted for in the boundary condition for 
the normal component of bulk velocity (|3|) and in the 
surface mass balance equation (j4|). The parameter r is 
the surface tension relaxation time. 

One would expect a generalised set of boundary con- 
ditions to have the no-slip condition as its limiting case. 
For the interface formation model this limiting case fol- 
lows from the limits p/(/3L), Ut/L — * 0, where L and U 
are characteristic length and velocity scales of the flow, 
applied to ©-([H]). Estimates for the phenomenological 
constants (3 and r, obtained from experiments on dy- 
namic wetting f25|, suggest that, say, for L ~ 10~ 2 cm, 
U ~ 1 cm s" 1 and p ~ 10 g cm -1 s _1 one has p/((3L) ~ 
10~ 5 and Ut/L ~ 1CT 6 . Then for flows which, unlike 
dynamic wetting, are not associated with infinitesimal 
length scales, boundary conditions ([2]H© to leading or- 
der in p/(/3L) and Ut/L, reduce to no-slip. 

Conditions ([2])-((6]) were derived using methods of irre- 
versible thermodynamics assuming that the equilibrium 
surface density p s el determined by the wettability of the 
solid substrate, is a constant, i.e. for a solid surface of 
a uniform wettability. This constant is related to the 
equilibrium contact angle that a liquid-gas free surface 
would form with the solid via the Young equation in the 
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following way. If we assume that the equilibrium surface 
tension of the solid-gas interface is negligible, the Young 
equation takes the form 

<?si = -crig cos6>, (7) 

where <j s i and <Ji g are the equilibrium surface tensions 
of the solid-liquid and liquid-gas interfaces and 9 is the 
equilibrium contact angle that the free surface forms with 
the solid. Thus, since according to ([6]), a s i = 7(P( ) — Pe) 
for a given liquid the Young equation J7]) allows one to 
express pi in terms of 9 (see (|Ti?)) below). 

In order to modify (J!])-© for a solid of variable wet- 
tability we must replace, in @ and ([5]) 

W^Vo- + p s F s , where F s — - — — . (8) 

Pi 

The reaction force F s acts on the liquid-solid intcrfacial 
layer from the solid surface and by balancing gradients of 
the equilibrium surface tension a e = a (pi), ensures the 
existence of a state of equilibrium. 

Consider the equilibrium surface density of the form 

Pi = 5 (PL + P & 2e) + 3 (P2e ~ Pie) ^ (*/0 , (9) 

where p\ e and p s 2e are the equilibrium surface densities 
in the far field as x — > — oo and x — > oo, respectively; i is 
the length of the transition region; a; is a Cartesian coor- 
dinate in the plane of the solid surface. The equilibrium 
surface density in the form of ^ allows one to investigate 
the role of a transition region between solids of different 
wettabilities and, after taking the limit l/L — > 0, find 
what boundary conditions for the surface variables one 
should use if this region is modelled as a solid-solid- liquid 
contact line. 

Finally, we will assume that the flow is plane-parallel 
in the (x, y)-plane of a Cartesian coordinate system, the 
origin of which is at the centre of the transition region, 
and that it is generated by constant shear of magnitude 
S in the far field, 

u — > S ij ■ r as r — > oo, (10) 

where i and j are unit vectors in the x and y directions 
and r is the radius-vector. Equations (fljl— pUj) now com- 
pletely specify the problem. 

It is important to emphasise that in the derivation of 
the model it is assumed that on the solid-facing side of 
the solid-liquid interface one has impermeability and no- 
slip. However, it is the velocity on the liquid-iacing side of 
the solid-liquid interface that is the boundary condition 
for the Navier-Stokes equations ([I]). In the classical no- 
slip condition it is assumed that there is no difference in 
velocity between the solid-facing and liquid-facing side of 
the interface, whereas in the interface formation model, 
the velocity on the liquid-facing side of the interface is 
determined by the interaction occurring in the surface 
phase and between the surface phase and the bulk. As 
a result, one can expect effective (or 'apparent') slip, i.e. 



the difference between the velocity on the liquid-facing 
side of the interface and the velocity of the solid surface, 
that is, in our case, u • (I — nn) ^ 0) and a non-zero 
normal component of velocity (a flux in/out of the surface 
phase, i.e. u • n / 0). 

It is convenient to non-dimcnsionalise equations 
(flOf using 

U = p~ 1 (Tl g , L = US~ 1 , P = pS, (Jig, P( 0) 

as the scales for velocities, length, pressure, surface ten- 
sion and surface density. Then in the bulk one has 

V-u = 0, Re (u • Vu) = -Vp + V 2 u. (11) 

whilst on the surface, where we use the notation u and v 
for tangential and normal components of velocity, 

( du dv\ \ ( da \p° dpl\ - 

\^j + d-x) + 2\Tx + ^llx-)=^ (12) 

v = Q (P s - Pi) , (13) 

^ --„-*>, («> 

a = X(l-p s ), (16) 

Pi = \ (Pie + P2e) + | [Pie ~ Pie) tanh . (17) 

and in the far field 

it— >l,u—>0 as x 2 + y 2 ^oo. (18) 

Here 

bp P O pTCTlg Gig 

(Tig (Jig P (0) 

It is noteworthy that, unlike the classical Navier condi- 
tion, in (fT2]l . which can be regarded as its generalisation, 
the tangential stress (the first term on the left-hand side) 
includes dv/dx. The deviation of p s from its equilib- 
rium value and the resulting adsorption/desorption lead 
to v 7^ via (fl"3"|) . Then, the spatial non- uniformity of 
this process makes dv/dx ^ thus 'switching on' this 
term in the tangential stress in (|12p and hence contribut- 
ing to the apparent slip. 
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The parameter Q in ([13)1 is the ratio of the charac- 
teristic mass flux into (out of) the liquid-solid interface 
associated with the adsorption (desorption) process trig- 
gered by the deviation of the surface density from its local 
equilibrium value and the characteristic mass flux due to 
convection in the bulk, whereas e _1 characterizes the ra- 
tio of the former and the characteristic divergence of the 
convective mass flux in the surface phase. In the present 
context, e is simply the product of the magnitude of shear 
in the far field S and the surface-tension-relaxation time 

T. 

Given that it is the equilibrium contact angle which 
a liquid-gas free surface would form with a solid, that we 
use as a measure of the wettability of a solid substrate, 
it is convenient to eliminate pf e , which are 'internal' pa- 
rameters of the model, in favour of di using ([7]) and (fTo]) : 



1 + A" 



COS ( 



(2=1,2) 



(19) 



Hereafter we will refer to the portion of the solid sub- 
strate with equilibrium contact angle 0\ and 62 as 'solid 
1' and 'solid 2', respectively. 

Following from [251 ] . where it was shown by analysing 
experiments on dynamic wetting that a ~ /3~ 1 , we will 
further assume that a/3 — 1 and hence that, in terms 
of our non-dimensional parameters, a/3 = 1. The analy- 
sis of experiments in [25j also provides estimates for the 
magnitude of phenomenological constants in the inter- 
face formation model's equations. Using these estimates 
and taking p ~ 1 g cm -3 , p ~ 1CP 1 — 10 2 g cm -1 s -1 , 
(Jig - 10 - 10 2 dyn cm" 1 , S ~ 10 3 - 10 5 s" 1 and 
I ~ 1CP 6 — 10~ 5 cm one arrives at a typical range of 
values for the magnitudes of the non-dimensional groups 



1(T 9 - 10 



Re 
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10~ 6 - 10~ 2 , 



e 

r~ io~ 6 

(* = 1,2) 
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As mentioned earlier, it is the double limit /3 _1 — * 0, e — > 
applied to equations (fP2")) - (fll)|) that results in the no-slip 
condition and hence all effects associated with deviation 
from the classical no-slip are at leading order in these 
parameters. 



III. SOLUTION 

The problem was solved numerically using the finite 
element method. Fig. [T] shows the streamlines of the flow 
for the case where solid 1 is more hydrophillic than solid 
2; the values of the dimensionless constants are given 
in the figure caption. As one can see, when the outer 
flow drives the fluid from a hydrophillic to a hydropho- 
bic zone, there appears a normal flux from the surface 
phase into the bulk. It is noteworthy that, as shown in 
Fig. [TJ the vertical component of bulk velocity is nonzero 
at y = 0, whereas for the classical Navier condition with 
different coefficients of sliding friction one has v = at 



y = [l2| and the normal component of velocity away 
from the solid appears solely due to the disturbance of 
the tangential flow at y = 0. 
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FIG. 1: Streamlines of flow over a chemically patterned sur- 
face in which values of the streamfunction i/i are given. Pa- 
rameter values Re = 0.01, e = 0.01, /? = 10 , Q = 100, F= 0.1, 
A = 10, 0i = 10° and 62 = 100° . 

In Figs. [2] [6] the distributions of the bulk velocity com- 
ponents along the solid surface are shown for the case in 
which solid 1 is more hydrophillic than solid 2. In each 
graph one can see (i) a positive normal velocity and (ii) 
variation in tangential velocity on the surface. The origin 
of the first of these effects is relatively straightforward. 
When fluid particles forming the interface are driven by 
the outer flow toward the region of lower equilibrium sur- 
face density, one has that in the disturbed equilibrium 
p s > pi, and hence, according to (fl3|) . v > 0. 

Importantly, it can be seen in Fig.[T]that the flux out of 
the surface phase occurs both in the hydrophillic (x < 0) 
and hydrophobic (x > 0) regions of the solid and extends 
itself well outside the transition zone. When solid 1 is 
more hydrophobic than solid 2 one observes the reverse 
effect, with the normal component of velocity directed 
towards the surface, corresponding to a flux into the sur- 
face phase. Once again this occurs on both sides of the 
transition region. 

In order to understand the distribution of the tangen- 
tial component of the bulk velocity on the surface we 
consider the terms on the left-hand side of the gener- 
alised Navier condition (fT2"[) which determine u{x, 0). For 
a plane-parallel shear flow over a homogeneous surface, 
i.e. far away from the wettability transition region, there 
is no deviation of the surface density from its equilibrium 
value (p s = pi) and hence no desorption from (or adsorp- 
tion into) the liquid-solid interface (v — 0). This is an 
obvious solution of (fTl"T) - (fT5|) . In dimensionless variables, 
the tangential stress on the surface, du/dy + dv/dx, has 
magnitude 1 (du/dy = 1, dv/dx = 0). 
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FIG. 2: Variation of the 'generators' of slip (the terms 
on the left hand side of the generalised Navier condition 
(|12[) ~) in response to a variation in the solid surface wetta- 
bility. In curve 1, T—[du/dy + dv/dx — 1] whilst in curve 
2, T=[da/dx + (\p s /p°) dp°/dx] /2. Results are obtained for 
parameter values Re — 0.01, e = 0.01, /3 = 100, Q = 1, 
F=0.1, A = 20, 0i = 10°, a = 100°. 



As one approaches the transition region, the situation 
changes. Curve 1 in Fig. [2] shows the distribution of the 
deviation of the tangential stress (the first term on the 
left-hand side in l|12p) from the value 1 across this region. 
Importantly, there are now contributions to the tangen- 
tial stress from both du/dy and dv/dx since the latter 
becomes nonzero due to variation in the normal velocity 
along the surface as a result of the spatially nonuniform 
desorption (adsorption) process caused by the deviation 
of p s from its local equilibrium value. As discussed ear- 
lier, this effect could follow neither from the standard 
Navier condition nor from any of its generalizations if the 
interface formation process and the associated mass ex- 
change between the interface and the bulk are not taken 
into account. 

The other factors that, according to (fT!?)) , lead to slip 
on the liquid-solid interface behave as follows. As men- 
tioned in the previous section, in equilibrium the surface 
force, [(Ap s /Pe) dp s e /dx\ /2, and the surface tension gradi- 
ent, (da/dx) /2, balance each other to ensure that there 
is no perpetual motion. However, when a shear flow in 
the far field is imposed, it disturbs the surface phase, and 
the above two terms are no longer in balance since they 
obviously depend differently on the variable surface den- 
sity p s . The deviation of the second term on the left-hand 
side of (TT2")) from zero is shown as curve 2 in Fig. [5] One 
can see that it is the variation in this second term in (|12|) 
that dominates and hence it is the imbalance between the 
surface tension gradient and the tangential surface force 
caused by shear flow that is mainly responsible for a vari- 
ation in slip on the surface as its wettability changes. The 
corresponding right-hand side of (|12p and the normal ve- 
locity for these parameter values are shown by curve 3 of 
Fig.H 



Fig. [3] illustrates the effect that the size of the param- 
eters e and (3 have on the solution. As was mentioned, 
all effects associated with deviations from no-slip, most 
clearly the nonzero normal velocity, are proportional to 
e and /3" 1 , and the results shown in the figure support 
this conclusion. Noticeably, the normal velocity on the 
surface appears to be symmetric about the centre of the 
transition region, where it achieves its maximum value. 
The tangential velocity seems to be antisymmetric about 
the origin, with a decrease in slip on solid 1 and increase 
in solid 2. 

Fig. d] shows the effect on the solution of the param- 
eters Re, Q and A. Curves 2 and 3, which are almost 
indistinguishable, demonstrate that the Reynolds num- 
ber has almost no effect on the dynamics of the surface. 
The local Reynolds number near the surface is always 
small due to the smallness of velocities in the vicinity of 
the solid surface even when the value associated with the 
global flow is relatively large, and one should indeed ex- 
pect inertial effects to have little impact on dynamics in 
the surface phase. 

The parameter Q determines the degree of mass ex- 
change between bulk and surface phase fT3")) . Curves 2 
and 4 in Fig. [4] show that the greater its value, the more 
pronounced the effect will be. By varying A one can see 
that the effect increases with the compressibility of the 
fluid. 

Fig. [5] shows that for a smaller transition region one 
has a sharper effect. In all cases the disturbance, i.e. 
the non-zero normal component of the bulk velocity and 
variation in slip, runs well outside the region of varying 
wettability. The results suggest that if instead of a re- 
gion, one treated the transition in wettability as a solid- 
solid-liquid contact line, then the disturbance to the bulk 
flow would be sharper than the effect associated with a 
finite transition region but would still occur over a finite 
region of the interface. 

It is interesting to note that although the shape of the 
curves in Fig. [5] differ, the integral of the normal velocity, 
i.e. the total flux out of the surface phase per unit time 

/oo 
v dx, (20) 
-oo 

which we consider as a measure of the effect that a 
patterned surface has on an adjacent flow, remains un- 
changed. The value of J is a sensible choice of measure as 
we have seen that it is the normal component of velocity 
that causes the noticeable deviation from plane-parallel 
shear flow. 

Given that the change in solid does indeed alter the 
flow of an adjacent liquid, consider how the magnitude 
of the effect is dependent on the choice of solids. This can 
be illustrated by looking at three different solids charac- 
terised by contact angles of 10°, 60° and 110° that a free 
surface would form with them. Fig. [6] shows the veloc- 
ity components for three different combinations of solids 
for the case where the first solid is more hydrophillic. 
The results suggest that the normal flux per unit time 
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FIG. 3: Distribution of u(x, 0) and v(x,0) for different e and 
fi. Curves 1-5 correspond to e = 0.01, = 500; e = 0.01, 
fi = 200; e = 0.01, fi = 100; e = 0.02, /3 = 100; e = 0.05, 
fi — 100, respectively; for all curves Re = 0.01, Q = 1, Z = 0.1, 
A = 20, 8i = 10° and 2 = 100°. 



is proportional to the difference cos 0\ — cos 82 which ex- 
plains the ordering in Fig. [BJ where we can see the case 
10° -> 60° gives a far smaller effect than 60° -> 110° de- 
spite the difference in contact angles being the same. The 
numerical analysis of the problem made it possible to ad- 
vance and then verify the following approximate formula 
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FIG. 4: Dependence of u(x, 0) and v(x, 0) on Re, Q and A. 
Curve 2 corresponds to Re = 0.01, Q = 1, A = 20; 1: A = 50; 
3 (indistinguishable within the graphical accuracy): Re — 10; 
4: Q = 2. Parameters e = 0.01, /3 = 100, f= 0.1, 6*1 = 10°, 
82 = 100° for all curves. 



for J: 



J 



eQ 

2f3X 



(cos — COS 62) ■ 



(21) 



This formula, which very accurately represent the nu- 
merical data over a wide range of parameter values, has 
been obtained as follows. The proportionality of the flux 
between the surface and bulk phases to the cosines of 
the contact angle is in fact what one should expect given 
that according to © and ([7]), cos 9 is a linear function 
of the equilibrium surface density, and it is the difference 
Pie ~ Pie that determines the flux. This conjecture of 
the linear dependence of J on cos Q\ — cos 02 is the only 
nontrivial step in obtaining (|21|) and it has to be veri- 
fied numerically. We also have to find the coefficient of 
proportionality in this dependence in terms of the dimen- 
sionless parameters of the problem. In order to do this, 
we used the following procedure. For the fixed values of 
#1 and 62, we chose a base state in terms of the remain- 
ing parameters and then varied one of these parameters 
away from its base state to consider the effect this vari- 
ation produces on J. After repeating this operation for 
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FIG. 5: Velocity on the surface for different widths of tran- 
sition region. Curves 1-3 correspond to I = 1, I = 0.1 and 
I = 0.01, respectively; for all curves Re = 0.01, e = 0.01, 
/3 = 100, Q = 1, A = 20, 0i = 10° and 2 = 100°. 



all parameters we arrive at equation (121|) which now has 
to be verified by independently varying all dimension- 
less parameters, including 9\ and 02. The accuracy with 
which this equation represents J is illustrated in Fig. [7] 

Order of magnitude arguments for the phenomcnolog- 
ical parameters and the analysis of experiments on dy- 
namic wetting [25f suggest that /3 ~ n/h, where h is the 
thickness of the interfacial layer (modelled here as an 'in- 
terface' of zero thickness). Using this estimate, for the 
dimensional flux per unit time out of the surface phase 
in a liquid/solid/solid system one has 
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FIG. 6: Velocity on the surface for different solid-solid pairs. 
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(22) 



Given that h is typically very small (a few nm for sim- 
ple fluids |26|), the above estimate highlights the subtle 
nature of the effects that we have described. 



FIG. 7: Dependence of flux on parameters Q, e, j3 and A. 
Parameters are varied around a base state, the results of which 
are represented by vertical crosses, of: (3 = 100, Q = 1, e = 
0.01, A = 20. Then, diamonds: e = 0.02; squares: Q = 2; 
triangles: (3 = 200; diagonal crosses: A = 40. For all curves 
0i = 90° whilst 02 is varied, and Re = 0.01 . Curves 1 — 3 
represent the predicted flux given by (21). 
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IV. CONCLUSION 

As was shown, the interface formation model applied 
to the flow over a solid surface of variable wettability 
allows one to describe the main features of this flow ob- 
served in molecular dynamics simulations, most notably 
the nonzero component of the bulk velocity normal to 
the solid surface. A natural link between 'wettability' 
interpreted in terms of the concept of the 'contact angle' 
featuring in the spreading of liquids on solid surfaces and 
a viscous flow over chemically-patterned solids with no 
free surface present has been established. Importantly, 
the interface formation model deals with these phenom- 
ena entirely within the approach of continuum mechanics 
with no artificial inclusion of intcrmolccular forces in its 



framework. An interesting feature that follows from the 
results is that slip, i.e. the difference between the tan- 
gential component of the fluid's velocity and the corre- 
sponding component of the velocity of the solid surface, 
results primarily from the disturbance of the force bal- 
ance in the 'surface phase' and not from the tangential 
stress, as follows from the standard Navier condition. 
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